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Abstract 



We analyze the finite temperature deconfining phase transition in 2+1 dimensional Georgi- 
Glashow model. We show explicitly that the transition is due to the restoration of the 
magnetic Z 2 symmetry and that it is in the Ising universality class. We find that neglecting 
effects of the charged W bosons leads to incorrect predictions for the value of the critical 
temperature and the universality class of the transition, as well as for various correlation 
functions in the high temperature phase. We derive the effective action for the Polyakov 
loop in the high temperature phase and calculate the correlation functions of magnetic 
vortex operators. 
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1 Introduction 



The deconfining phase transition in QCD has been a subject of numerous studies in recent 
years. Some aspects of the physics of the high temperature phase appear to be perturba- 
tive and can be studied in a controlled way at asymptotically high temperatures. However 
in the phase transition region itself the QCD coupling is large and the physics is domi- 
nated by the nonperturbative soft sector. Analytic understanding of this region is therefore 
extremely difficult. It thus seems useful to find a simpler toy model where a similar decon- 
fining phase transition is within the weak coupling regime and can be studied analytically. 
Optimistically one can hope to learn some lessons from such a model which could be of 
use in our understanding of real life QCD. 

A model that is very suitable for this purpose is the Georgi-Glashow model in 2+1 
dimensions - Yang Mills theory with a Higgs field in the adjoint representation. Over twenty 
years ago Polyakov [|IJ showed that in this theory the instantons (monopoles) dominate the 
nonperturbative dynamics and disorder the Higgs vacuum leading to a linear confinement 
of electric charges even at arbitrarily weak coupling. This was the first demonstration of 
the dual Meissner effect (or dual superconductivity) in a non-Abelian gauge theory where 
monopoles condense and expel electric fields from the vacuum, providing an example of 
the enigmatic confinement phenomenon in controllable circumstances. At zero temperature 
the 3D Yang-Mills-Higgs system is therefore always in the confining phase, even though 
classically the confining and Higgs regimes appear to be separated by a phase transition. 

In this paper we will study the theory at finite temperature with particular emphasis 
on the deconfining phase transition. We note that a similar study appeared in the literature 
several years ago 0. Our results are both qualitatively and quantitatively different from 
those of 0. In particular we show that the deconfining phase transition is primarily 
due to the plasma effects of charged excitations (W bosons), an effect that was neglected 
completely in the analysis of 0]. As a result our predictions for both the value of the 
critical temperature and the universality class of the phase transition are different. 

The structure of this paper is the following. In section 2 we review the model and 
discuss the two complementary views of the confinement mechanism: one based on the 
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monopole plasma |TJ and the other based on vortex condensation |3|, |j. In section 3 we 
discuss the naive expectations for the deconfmement phase transition which follow from 
these two approaches. We show that these expectations are different and explain why the 
monopole plasma intuition (the basis of the approach of 0) should be refined in order to 
discuss the phase transition properly. In section 4 we describe the renormalization group 
analysis of the phase transition. We determine the value of the critical temperature, the 
universality class and discuss the crucial role played by the electrically charged particles. 
In section 5 we relate our analysis to the standard tool used in the high temperature phase 
- the effective potential for the Polyakov line. We derive this effective potential in the 
deconfining phase and show with its help explicitly that the correlation function of vortex 
operators at T > T c decreases exponentially at large distances and thus that the magnetic 
Z 2 symmetry is restored. Finally in section 6 we discuss briefly our results. 

2 The model 

We consider the £77(2) gauge theory with a scalar field in the adjoint representation in 
2+1 dimensions. 



£ = J d 3 xtr (F IM ,F fJ,v ) + J d 3 x 



±(D„h a ) 2 + ~(h a h a -v 2 ) 2 



(1) 
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We adopt the notation = ^A a ^T a , F^ = d^A u - d v A^ + [A^A^], h = \h a T a , and 
Dfji = dfji + [Ap, h] and the trace is normalized as tr(r a r b ) = 2S ab . 

We will be interested throughout this paper in the weakly coupled regime v ^> g 2 . 
In this regime perturbatively the gauge group is broken to £7(1) by the large expectation 
value of the Higgs field. The photon associated with the unbroken subgroup is massless 
whereas the Higgs and the other two gauge bosons are heavy with the masses 

M 2 H = 2Xv 2 , Ml = g 2 v 2 . (2) 

Thus perturbatively the theory behaves very much like electrodynamics with spin one 
charged matter. The nonperturbative effects however are very important at large distances. 
Their effect is that the photon acquires finite mass and the charged W ± become linearly 
confined at large distances with nonperturbatively small string tension. 



2.1 The monopole plasma 

The nonperturbative configurations that dominate the generating functional in the weakly 
coupled regime are solutions to the classical Euclidean equations of motion, 

= [h, D v h\ , 

D^D^h = -\(h a h a - v 2 )h . (3) 

The perturbative vacuum manifold is SO(3)/U(l) and this has a nonzero second homotopy 
group, U 2 (SO(3)/U(l)) = Z. The theory thus allows topologically non-trivial solutions. 
Polyakov showed [|l[] that, in the vacuum only far separated monopole configurations are 
relevant. The unit charge configuration is of the well known 't Ho oft- Polyakov monopole 
form, 

h a (x) = x a h(r) 

A ffi = I l^ux u (l - 00 + 8 a y 2 + (rA - 2 )x%] (4) 
where x a = x a /r. The non-Abelian field strength is 

f; v = ^ fe ^(</>? + ^ - 1) + ^ 

+(5 au x^ -5 a ^x u )-(<f)' -AM, (5) 
r 

where prime denotes differentiation with respect to 3D radial coordinate r. The finiteness 
of the action (in the gauge A = <p 2 = 0) requires h(0) = , h(oo) = v, 4>i(0) = 1, and 

<M°o) = o. 

The monopole has a core size of the order of the vector meson mass. Outside the 
monopole core, the theory is Abelian and one can define the Abelian field strength as, 

F, v = ^F; u - ±e abc h a (D^h) b (D„hy. (6) 

The monopoles (instantons) interact via a Coulomb potential in the Euclidean three 
dimensional sense, and the partition function of the Coulomb gas can be written as a path 
integral over the "dual photon field" <fi. The effective Lagrangian for the field <fi reads 



The dynamically generated photon mass, due to the Debye screening of the monopoles, is 
given in terms of the fugacity of the monopoles as 



M = j^. £ = constant— 1 
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e(^g-) is such that 1 < e < 1.787 §, and e(oo) = 1. The field <p is by far the lightest 
excitation in the theory and thus eq.(|7|) is the low energy effective Lagrangian of the theory 
valid below the energy scale M\y. The fundamental Wilson loop calculated in this effective 
theory |l| has an area law fall off, indicating confinement. 

2.2 The vortex condensation 

A complementary view of confinement in this theory is the spontaneous breaking of the 
magnetic Zi symmetry due to the condensation of magnetic vortices 0. 

Ignoring monopole effects, the model has the magnetic U(l) symmetry generated by 
the Abelian magnetic field. The conserved current is the dual field strength 

F» = \^ v xF uX . (9) 

As discussed in detail in [|J], due to the monopole contributions this symmetry is anomalous, 
and only the Z% subgroup is conserved. The non-anomalous Z2 magnetic symmetry group 
is generated || by the large spatial Wilson loop which encloses the system 

W(C -> 00) = exp^ J d 2 xB(x) (10) 

The order parameter for this Z2 transformation is the operator that creates an elementary 
magnetic vortex of flux 2n/g [||, ||] 

2iri r h a 

K( x )=exp— j^—E^x). (11) 

This operator is local, gauge invariant and Lorentz scalar. Together with the spatial Wilson 
loop it forms the order-disorder algebra 

W(C -> oo)V(x) = -V(x)W(C -> 00) (12) 
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It can be shown || that the vortex operator condenses in the vacuum and its expecta- 
tion value is determined by the gauge coupling constant < V > 2 = g 2 /8ir 2 . The low energy 
theory from this point of view is given by the Z2 invariant theory of the vortex field V 



s\ 2 , M 2 
lix 2 > 4 



£eff = d»V d»V* - X(V V* - f^) 2 + — {V 2 + (V*) 2 }. (13) 



The vortex selfcoupling in this effective Lagrangian can be determined as || 

A = 5-^. 14 

9 

Note that as a low energy Lagrangian, eq.(|T3"D is indeed consistent with eq.(^). At 
weak gauge coupling the quartic coupling A is very large. In this nonlinear a-model limit 
the radius of the field V is therefore frozen to its expectation value. The only relevant 
degree of freedom is the phase 

V(x) = -jL-e* (15) 
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Substituting this into eq.(|13"D one indeed obtains precisely eq.(|7D with the identification 
= 2x- It was shown in || that this identification is indeed the correct one, and the factor 
2 follows from the fact that the magnetic charge of the monopole is twice the fundamental 
value relative to the electric charge of the elementary charged excitation W : gu9e = 47r. 

The main difference between eq.(|13D and eq.(|7]) is in the description of the charged 
sector. As explained in ||, the charged particles are solitons of the field V with unit 
winding number. In the nonlinear cr-model limit those are just vortices of the phase X- 
Polyakov's effective Lagrangian eq.([7p does not allow such vortex configurations. In the 
vortex configuration the field x (and <p) is discontinuous along an infinite cut. The cut 
contributes to the kinetic energy term in eq.(|7]), and this contribution is both ultraviolet 
and infrared divergent. The ultraviolet divergence is of no importance by itself, since 
the theory is defined with the intrinsic UV cutoff Mw- However the infrared divergence 
indisputably puts description of the charged states beyond eq.(|7|). On the other hand 
eq . (|i~3"D recognizes the fact that the field x is a phase, and therefore 27r discontinuities in it 
are unphysical and do not cost energy. Thus vortex configurations in eq.([T3|) are allowed. 
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Although they are confined dynamically with a linear potential [|J, the string tension for 
this potential is a dynamical effect proportional to £}l 2 and has nothing to do with the 
extraneous cut contribution. 

In fact, if anything the Lagrangian eq. ([l3|) underestimates the energy of such a vortex. 
The core energy of the vortex should be equal to the mass of the charged vector boson 
My/, since this is the lightest charged excitation in the theory. On the other hand the core 
mass of a vortex in the Lagrangian eq. ([13]) is given by the UV contribution of the Coulomb 
potential. With the cutoff of order My/ this is of order g 2 In My/ / g 2 . This is not surprising, 
since the mass of the W boson indeed comes from the distances of order of its Compton 
wave length, and is thus well inside the UV cutoff of the effective theory. The situation can 
be improved by adding to the effective Lagrangian a higher derivative term of the Skyrme 
type 

5L = A(£^ u \d u V* d\V) 2 (16) 

with 

9 My/ 

With this extra term the core energy of the vortex indeed becomes of the order of My/. 

As we will see in the following, the correct description of the charged sector is crucial 
for the understanding of the phase transition^. 

3 Deconfining phase transition - The Rough Guide 

In this section we give a naive discussion of the expected nature of the deconfining phase 
transition based on the monopole versus vortex description of the low energy theory. 



5 We note that this is not the only instance in which the charged states play important role. For example 
as discussed in H the presence of these states leads to the breaking of the string of charge two, which in 



the Polyakov's effective Lagrnagian eq.)|^ is strictly stable. 
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3.1 The monopole binding 

Since at zero temperature the monopole contributions are the only relevant ones for con- 
finement, one may be tempted to assume that also at finite, low enough temperature all 
other effects are not important. The expected value of the transition temperature, as we 
shall see in a moment is of order g 2 , which is much smaller than any mass scale in the 
theory except for the photon mass. One could therefore start with the working hypothesis 
that the single self-interacting photon field (or equivalently the monopole ensemble) should 
be a valid description of the phase transition region. This is the point of view taken in , 
and we now briefly describe the conclusions it leads to. 

The first thing to note is that at finite temperature the interaction between monopoles 
is logarithmic at large distances. The reason is that the finite temperature path integral 
is formulated with periodic boundary conditions in the Euclidean temperature direction. 
The field lines are therefore prevented from crossing the boundary in this direction. The 
magnetic field lines emanating from the monopole have to bend close to the boundary and 
go parallel to it. So effectively the whole magnetic flux is squeezed into two dimensions. 
Qualitatively the situation is shown in figure 1. 

A 



Figure 1: The field of a Monopole- Instanton at finite temperature. 
The length of the time direction is (3 = 1/T, and thus clearly the field profile is two 
dimensional on distance scales larger than f3. Two monopoles separated by a distance larger 
than f3 therefore interact via a two dimensional rather than a three dimensional Coulomb 
potential, which is logarithmic at large distances. Since the density of the monopoles is 
tiny p M oc £, already at extremely low temperatures T oc £ 1//2 the monopole gas becomes 
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two dimensional. The strength of the logarithmic interaction is easily calculated. The 
magnetic flux of the monopole should spread evenly in the compact direction once we are 
far enough from the monopole core The field strength should have only components parallel 
to the spatial directions. Since the total flux of the monopole is 2ir/g, the field strength 



far from the core is Fj = ^^j, and thus the strength of the infrared logarithmic interaction 

-*2 I „1 



is T jg 

It is well known [0] that the two dimensional Coulomb gas undergoes the BKT phase 
transition. In the usual Coulomb gas, where the strength of the interaction A does not 
depend on temperature, the particles are bound in pairs at low temperature and unbind 
at high temperature T > Tbkt = 2tt\, where the entropy overcomes the energy. This 
is the standard BKT phase transition |7[ which determines the universality class of U(l) 
symmetry restoration in 2 dimensions. In the present case the situation is reversed, since 
the strength of the interaction itself depends on the temperature. At low temperature the 
interaction is weak, and therefore the particles (monopoles) are free. As the temperature 

T 2 

grows, the interaction becomes stronger until at Tbkt = 27r ^f T the energy overcomes the 
entropy. Above this temperature the monopoles bind into neutral pairs. Thus based on 
this simple picture we expect the theory to undergo a BKT phase transition at 

g 2 

Tbkt = tt~ 0-°) 
in 

Below this temperature the photon should be massive, while above this temperature it 
should be massless since the cosine term in the Lagrangian eq.(|7|) is irrelevant. 

This is precisely the logic followed in M. One can go further and perform more 
quantitative calculations using the dimensionally reduced Polyakov effective Lagrangian. 
Dimensional reduction should be perfectly valid in this theory since the critical temperature 
is much larger than the photon mass. Thus one can consider the two dimensional sine- 
Gordon theory 

c=£rMr + (19) 

This theory, in full agreement with the previous discussion has a BKT phase transition at 
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the temperature given by eq.([TJ), above which the "dual photon" field becomes massless. 

Moreover, since the 2-d sine-Gordon theory is exactly solvable, one can calculate 
various correlation functions. In particular the authors of || calculated the string tension, 
and found that it exhibits quite strange behavior. At low temperatures it follows an 
expected pattern - that is, it decreases with temperature. It becomes extremely small at 
temperature much smaller than Tbkt] however just before Tbkt it starts rising sharply 
and at T BKT actually diverges. 

Although the logic leading to the previous discussion seems quite sound, there are 
several puzzles that arise regarding the results. 

• The photon becomes massless in the high temperature phase. In other words, the 
correlation length in some physical (gauge invariant) channel becomes infinite at high 
temperature. This contradicts our physical intuition, since one expects that at high tem- 
perature any physical system described by a local field theory should become maximally 
disordered with zero correlation length. 

• The phase transition is of BKT type, and therefore belongs to the U(l) universality 
class. On the other hand we know that the global symmetry that is restored at the phase 
transition is the magnetic Z 2 f|, and we expect the universality class to be 2D Ising. 

• The divergence of the string tension at the critical point looks unphysical. 

• Finally, as we explained in the previous section, the effective theory eq. (|7D does not 
allow charged states. It is then difficult to understand in what sense the high temperature 
phase can be viewed as deconfined. 

3.2 Magnetic symmetry restoration and the charged plasma 

Let us now take eq. fll3|) as our starting point. What does one expect from the phase 
transition in a simple scalar theory of this type? Firstly, clearly one expects the criti- 
cal temperature to be of order of the expectation value of the scalar field, and therefore 
parametrically of order g 2 . At finite temperature one expects the generation of a positive 
thermal mass proportional to the temperature and to the coupling constant. Thus the 
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thermal contribution to the effective potential should be of the form 



5 T L = x\TV*V (20) 

At T = g 2 /4ir 2 x the total mass term becomes positive, the VEV of V vanishes and the phase 
transition occurs. One can estimate this temperature more precisely using the following 
simple argument. Let us neglect for now the monopole induced term. Then we are dealing 
with the XY model at finite temperature. The dimensional reduction should again be a 
valid approximation, and thus essentially we have to analyze a 2D model. Now a 2D XY 
model can be mapped into a sine-Gordon theory of a dual field x- Performing this dual 
transformation^ , we find the Lagrangian 

L = ^-(d l x) 2 + f icosx (21) 

where \x is the fugacity of the vortices in the XY model. This sine-Gordon theory has a 
phase transition at 

Txy = f- (22) 

Therefore we may expect the magnetic symmetry restoring phase transition at Txy- Al- 
though parametrically this temperature is of the same order as the BKT transition tem- 
perature discussed in the previous subsection, it is four times lower. Another important 
difference is that the nature of the phase transition is in fact completely different. The 
phase transition in the model eq. (pip is due to the unbinding of vortices of the field V. 
Above Txy the vortices are in the plasma phase. These vortices, are as discussed earlier 
precisely the charged W bosons of the original Georgi-Glashow model. Thus this phase 
transition is just what one would naturally call the deconfining phase transition and the 
vacuum above Txy is the charged plasma. 

Note that in this discussion we have neglected completely the effect of monopoles, that 



is the last term in the Lagrangian eq. (|13|). Interestingly enough we see that the plasma 

6 To fix the normalization of the kinetic term we should bear in mind that a vortex corresponds to a 2ir 
jump in the field \ — 2</> rather than the field (j> in eq. ([?]). See also the discussion following eq. (|J 
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phase is reached at a temperature which is much lower than the Tbkt discussed earlier, 
and thus the monopole binding is irrelevant for the dynamics of the deconfmement. This is 
not to say that the presence of the monopoles is irrelevant altogether. Clearly, omitting the 
monopole induced term we enlarged the symmetry of the system from Z 2 to U(l). Hence 



the effective XY model description and the U(l) universality class predicted by eq.(21). 
The analysis as it stands is correct for noncompact electrodynamics with charged matter, 
but not for the Georgi-Glashow model. Next we remedy this problem. 

The effect of the monopoles can be qualitatively understood in a simple way. Their 
presence leads to a confining potential between the charges, that is a linear interaction 
between the XY model vortices. Thus the monopoles suppress the variation of the phase 
of the field V except inside the confining string, the width of which is given by the inverse 
photon mass d = 1/M oc £ 1/2 . Whether the presence of the monopoles is important at 
the "would be" phase transition point Txy depends crucially on the density of charged 
particles. If the density of the charged particles is very high, so that the average distance 
between them is smaller than d, the presence of the monopoles is immaterial since the 
phase of V is disordered already on the short distance scale. However if the density of 
charges is low, so that the distance between them is larger than d, the presence of the 
monopoles will suppress the phase transition. In this case at Txy the vacuum will not be 
disordered, but will look like a dilute gas of charged particles with strings between them. 
The actual phase transition will then occur at a higher temperature, where the average 
distance between the charges equals the inverse photon mass. 

Interestingly enough in the present case this temperature turns out to be twice the 
value of Txy and therefore still much lower than Tbkt- We will show this rigorously in 
the next section. Qualitatively though this is easy to understand. The density of charges 
is proportional to their fugacity 

H oc e t (23) 
This should be compared to the square of the photon mass, which is given by the zero 
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temperature monopole fugacity. In a theory with very heavy Higgs 



£ oc e~^~ (24) 

The two are equal (up to subleading corrections) at 

Too = £ (25) 

More generally, when the Higgs is not infinitely heavy, the monopole fugacity is smaller 
and thus the transition temperature will be lower 

T ^4&) (26) 

For lighter Higgs, T c gets closer to Txy, but is always greater so that Txy < T c < Tbkt- 

In the next section we will present a more complete analysis based on the renormal- 
ization group and exact bosonization. This analysis confirms the simple picture presented 
here. Thus all the puzzles raised in the previous subsection simply disappear. 

• The photon never becomes massless. Even without the monopoles in the plasma 
phase it acquires the Debye mass given by the cosine term in eq. (pi]). This mass rises 
with temperature, and thus the physical correlation length decreases. 

• Since the monopole term is still relevant at the phase transition, the universality 
class must be Z 2 - We will see this explicitly in the next section. 

• The analysis of the previous subsection is only valid below Txy, and thus the 
divergence of the string tension at T BKT > T X y has nothing to do with the physics of the 
Georgi Glashow model. 

• Finally, the phase transition is driven primarily by the unbinding of charged particles 
and thus indeed has all the flavor of a deconfinement transition. 

Let us now turn to a more quantitative analysis of what is stated above. 

4 Deconfining phase transition - renormalization 
group analysis 

For a more formal analysis of the phase transition we find it convenient to use the sine- 
Gordon formulation in terms of the phase field x modified to take explicitly into account 
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the finite probability of the appearance of vortices. This Lagrangian has the form 

C = -f— (d, X ) 2 + C cos 2 X + fi cos x (27) 
where ( is related to the monopole fugacity ( = £/T and X is the field dual to \i 

id,X = l^f^X- (28) 

4.1 The Lagrangian 

The way to derive this Lagrangian is the following. The partition function of the sine- 
Gordon model in the presence of one vortex is 

Z(x) = J D[ X }exp{-J ^y^id.x- My, x)) 2 + Ccos2 X } (29) 

The "external current" is 

Uy,x) = 2imM6{yeC) (30) 

with C a curve that starts at the location of the vortex (the point x), and goes to infinity, 
and n M is the unit normal to this curve. The insertion of this current forces the derivative 
of x to have a discontinuity across the curve C, so that x jumps by 2n. This forces 
X to have one unit of vorticity concentrated at the point x. Note that even though 
explicitly depends on the curve C, the partition function itself does not, since changing 
the integration variable 

X (x)^x(x)+2tt, xeS (31) 

where the boundary of S is C — C is equivalent to changing C into C in the definition of 
the current. The extra linear term in the exponential in (eqj29|) is 

f_ 

2-kT Jc 

which is equivalent to (eq.|28D. An antivortex at y obviously is created by — j M . To create 
several vortices one just inserts an external current which is the sum of the currents which 
create individual vortices. 
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IX 



/ dx^e^dyX- (32) 
J c 



A dilute ensemble of vortices and antivortices with (small) fugacity \i is then given by 



The summation over the number of vortices and antivortices can be easily performed 
leading to the partition function with the Lagrangian eq. (WA) - The constant fi is the 
vortex fugacity scaled by the effective UV cutoff imposed on the integration over the 
coordinates. The vortex fugacity of course is none other but the fugacity of the charged 



The cutoff a is related to the Compton wave length of the W boson, but a more careful 
determination of it should take into account the fact that in the process of dimensional 
reduction all modes with frequencies above T have been integrated out. We will not 
attempt the determination of a, but only note that it is some combination of the scales 
Myj/ and T oc g 2 , and as such its value always plays a role secondary to the exponential 
factor of fugacitj|]. 

An alternative way to derive the Lagrangian eq. ( p?|) is to start directly from the 
effective theory eq. (|T3| ) with the extra Skyrme term eq. (|T6D. In the nonlinear cr-model 
limit one can cleanly separate the phase of the field V into a smooth part x an d the vortex 
contribution. The Skyrme term then is proportional to the energy of one vortex, and with 
the Skyrme coupling chosen the way we discussed in Section 2 is just equal to Mw for a 
one vortex configuration. The dilute vortex gas approximation then reproduces again the 
Lagrangian eq. (W^i)- 

4.2 The renormalization group equations 

Our starting point for the discussion of this section is eq. (p7|). Since both £ and /i are 

small, the importance of different terms in the Lagrangian is determined by their respective 

7 For the same reason the effective value of the monopole induced mass £ also changes in the dimen- 
sionally reduced theory. This change has been calculated in Again we will not dwell on this finite 
"renormalization" since this effect is subleading in the same sense. 




(33) 



W, 



/i = a e t 



(34) 
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conformal dimensions calculated in the free theory. The total conformal dimensions of the 
operators : cos 1\ '■ , '■ cos % '■ respectively are, 

* = A„^. (35) 

Thus at low temperature the charges are irrelevant (A M > 2) and the monopoles are 
relevant (A^ < 2) and the theory reduces to the sine-Gordon model. At the temperature 
T = g 2 /4:7r the two operators become equally relevant since their conformal dimensions are 
equal 

\ = \ = 1- (36) 

The phase transition point can be determined from the structure of the fixed points 
of the renormalization group equations. The critical point is the IR unstable fixed point 
of the RG flow. The renormalization group equations of the sine-Gordon theory have been 
studied perturbatively in the literature both in the absence and in the presence of the 
vortices fTI|| . In terms of the dimensionless parameters 

47T 

t = —T, fi = fia 2 , ( = (a 2 (37) 

g 2 



the lowest order RG equations read 

dt 



dX 

The fixed point structure of these equations is simple. 
1. The point T 



2 = n 2 (f-t 2 C 2 ), (38) 
f = (39) 
= (2 - t) C (40) 



t = 0, /t = 0, C = oo (41) 



is clearly the zero temperature fixed point. Here the long distance physics is dominated by 
the monopole induced mass term. 

16 



2. The point T c 



oo 



t — OO, 



/I = OO 



C = o 



(42) 



is the high temperature fixed point. Here the infrared properties are detemined by the 
charged plasma effects. 



This is the IR unstable fixed point. This is precisely the critical point that corresponds to 
the deconfining phase transition. At this point both the monopole and the charge plasma 
induced terms are equally important. 

Note that the fixed line t > 2, ( = of massless theories which corresponds to the 
"confined monopole plasma" is not present, since for t > 2 the charged plasma induced 
mass term is strongly relevant. The same is true for the would be fixed line t < 1/2, jl = 0, 
which is present in the absence of the monopoles and describes the low temperature massless 
phase of noncompact electrodynamics. 

It is instructive to see how the RG equations formalize our qualitative arguments of 
the previous section. In particular they make clear the role of the point t — 1/2 which 
in the absence of monopoles would be the point where the charged induced term becomes 
relevant. If one starts the evolution from the initial condition t — 1/2, jl the running 
temperature will increase, and jl will grow into the infrared, while ( will grow for a while 
until the running temperature reaches two. From this point on ( will decrease and the 
system will approach the high temperature fixed point. This corresponds to the situation 
where at the "would be" critical point t — 1/2 the density of charged plasma is so high 
that the mean distance between the particles is smaller than the width of the confining 
string. 

On the other hand if the initial condition is jl < (, the temperature starts decreasing, 
making the coupling jl immediately irrelevant. Thus ji monotonically decreases to zero 
and practically does not affect the flow of the other couplings which steadily flow to the 



3. The point T GG 



t = l, 



jl = (, jl = oo 



(43) 
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zero temperature fixed point. This is the situation where at t — 1/2 the density of charged 
plasma is low, which is indeed true in our model. 

Interestingly enough, at t = 1 the initial conditions in our model in the case of very 
heavy Higgs are such that the RG flow starts almost exactly in the region of attraction of 
the fixed point Tgg- At this value of temperature the monopole and charge fugacities are 
equal. The only difference between ft and ( then is in the prefactors, which are possibly 
different combinations of My/ and g 2 . If ft and £ were exactly equal, the only thing that 
happens along the flow is that their values grow, but t does not change. Due to the small 
difference their values are equal at a temperature which slightly differs from 1 

^^ + 0(1 ), T-£n +0 (£,] (44, 

Thus the initial temperature for which the system is in the region of attraction of Tgg is 
slightly different. We conclude that the critical temperature of the Georgi-Glashow model 
is indeed, up to corrections of order g 2 /M w , given by T GG = 

A comment is in order about our use of the perturbative renormalization group equa- 
tions. As written the equations (^) are only valid for small ft and (. The initial values 
of these constants in our model are indeed small, and thus the flow is well described by 
these equations for quite a while. The fixed points however correspond to large values of at 
least some parameters, and thus quantitative details of the flow in this region of parameter 
space are different. However the existence of the three fixed points T , and Tqg is an 
exact statement. For T and this is a consequence of the fact that one of the couplings 
vanishes (/2 or Q and thus these fixed points are the same as in the corresponding sine- 
Gordon theory. The existence of Tgg is a consequence of the exact duality in the model eq. 
(p7|). The transformation 2\ — > X m the path integral leads to the same partition function 
but with new parameters as follows t — > 1/t, p, — > £ — ► ft. This symmetry is nothing 
but the famous T-duality (see for example [|12[ and references therein). The point Tgg is 
the fixed point of this transformation and is thus guaranteed also to be a fixed point of the 
exact renormalization group equations. 
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4.3 The universality class of the phase transition 

Since the deconfining phase transition is due to the restoration of Z2 symmetry we expect 
that it is in the 2-dimensional Ising universality class. This can be shown exactly by 
studying the theory eq. (p7|) at the fixed point t = 1. The following discussion follows 
closely jnj[See also chapter 21 in Ipfl l- 



The theory described by eq. (p7|) can be fermionized by using the standard bosoniza- 
tion/fermionization techniques. Since at t = 1 both cosine terms have dimension 1, the 
resulting fermionic theory is a theory of free massive fermions. In our notations, the Dirac 
fermionic field is defined as 

i/)r = a' 1/2 iexp[i(x + |)], ip L = a~ 1/2 exp[-z(x - |)], . (45) 



The kinetic term in eq. fl27|) then becomes the kinetic term of the field ip, while the cosine 
terms become 

aT x cos 2% = ityj^L — H.c] 

a~ x cosx = i[ip R ip L — H.c.]. (46) 

Thus the mass term in the fermionized Lagrangian is diagonalized by introducing the real 
Majorana fermions 

" = — • ° = ^J2- (47) 

The mass of the fermion p is fia + (a, while the mass of a is fia — (a. For \i = ( the model 
contains one massive and one massless fermion. At the fixed point jx —>■ 00 the massive 
fermion decouples. Thus at the point Tqg the theory is that of one massless Majorana 
fermion. It is well know that this theory precisely describes the critical point of the 2D 
Ising model. 

Thus we indeed see that the phase transition is in the universality class of the 2D 
Ising model. 
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5 From the magnetic vortices to the Polyakov line 

It is standard practice to study the high temperature phase of gauge field theories using 
the effective action in terms of the zeroth component of the vector potential Aq. In this 
section we want to relate our analysis which was performed for the phase of the vortex 
operator x to this standard procedure. 

The relation in fact is not difficult to make. The free energy of the charged bosons 



in the Lagrangian eq. (|2"7 ) is given by the term cosx- On the other hand this same free 
energy in the standard calculation is represented by the insertion of the Polyakov line with 
charge two. We thus identify the Polykov line with the exponential of the dual field 

P = exp{~x} (48) 

and the dual field \ with the Abelian vector potential 

X = 2gf3A (49) 

Up to the monopole induced term, the Lagrangian eq. ( p7[) is equivalent to the sine-Gordon 
theory of the field x i n ec L- (HP- With the identification in eq. (|9|), this is 

|(^o) 2 + /icos(^A ) (50) 

The monopole induced term also can be written in terms of the vector potential. Its form in 
terms of the Polyakov line P is precisely the same as the form of the charge plasma induced 
term in terms of the vortex field V. Recall that the origin of the plasma induced term 
cosx in eq. (]27D is the dimensionally reduced "Skyrme term" of eq. (|l6|) . Thus just like we 
derived the sine Gordon Lagrangian eq. fl27D starting from the effective Lagrangian for the 
vortex field, we can follow the same steps backwards but this time expressing everything in 
terms of P. Due to this duality between V and P we conclude that the monopole induced 
term is the dimensionally reduced Skyrme term for the Polyakov line. The full effective 
Lagrangian for the vector potential is then 

L = |(<9^o) 2 + /i cos(|Uo) + ^ ln (0 ^APd 3 P*) 2 (51) 
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The monopole in this Lagrangian is represented by a vortex of the field ^Aq with unit 
vorticity. The coefficient of the Skyrme term in eq. (]5lD is such that the action of such a 
unit vortex is equal to the action of the "core" of the monopole. 

The cosine term in this expression is the potential for A induced by the nonvanishing 
density of the charged particles in the thermal ensemble. Naturally, it contains the Debye 
"electric" mass term for Aq and also higher interactions. Note that as opposed to strongly 
interacting theories, where similar effective Lagrangians have been derived only in the 
derivative expansion, the Lagrangian eq. ([H]) is valid on all distance scales longer than 1/T. 
Thus in principle it can be used to calculate correlation functions in a large momentum 
range. However, if one is interested in the long distance behavior of the correlators, at 
temperatures above Tqg the higher derivative Skyrme term can be neglected. This is 
in accordance with our analysis of the previous section which showed that the monopole 
induced term is irrelevant above Tqq. 

An interesting example of a correlation function is the correlator of vortex operators 
< V(x)V*(y) >. At low temperatures, since the magnetic Z 2 symmetry is spontaneously 
broken this correlation function tends to a constant at large separations. The corrections 
to this constant value are given by the correlator of the field x- At zero temperature the 
correlator is 

q 2 1 q 2 16tt 2 

< V(x)V*(y) >= ^ exp{-- < x(x)x(v) >} = exp{ — ^-Y 3 (x - y)} (52) 

where Y^(x—y) is the 3d Yukawa potential with the mass M. At low temperature (T < Tgg) 
the infrared asymptotics (at distances \x — y\ > 1/T) is instead 

< V{x)V\y) >= £^eM-^-Y 2 (x-y)} (53) 

with Y 2 (x — y) - the 2D Yukawa potential with the mass which includes the effects of 
integrating out the nonzero M atsubara modes, as calculated in |2| . This expression follows 
from eq. (1271) neglecting the cos x term, which is indeed negligible in the infrared. 



At high temperatures T > T GG instead it is convenient to use eq. (]5lD with the 
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omission of the third termQ The calculation of the correlator of the vortex operators 
proceeds along the lines of |T5|, The insertion of V(x) and V*(y) creates the Z 2 domain 
wall stretching between the points x and y. In terms of the sine-Gordon theory eq. ( |51~D 
this domain wall is just the kink, and thus the domain wall tension is equal to the soliton 
mass M, [fill. The correlator then is 



< V(x)V*(y) >= exp{-M s |x - y\} 



(54) 



with the soliton mass reading as 



M. 



-i 2r(p/2) 
1 yftrm 



v p+i > 
rc-M 



2/2 



p+i 

2 



(55) 



where 



V 



r 



(56) 



8ttT - g 2 ' 

Thus we find that the correlation function of the vortex operators decreases exponen- 
tially in the high temperature phase as it should in the phase with restored symmetry. 
Note that were we to neglect the Debye mass term we would have found that below Tbkt 
the correlator tends to a constant at infinity, while at T > Tbkt it does decay at large 
distances but only as a power, since the mass of the soliton would vanish. 

We note that the result eq. ( |54"D implies that in the deconfining phase the monopoles 
are confined linearly. We remind the reader that insertion of an operator V is equivalent to 
creation of a monopolef]. Thus the logarithm of the correlation function can be interpreted 
directly as the intermonopole potential. Indeed at zero temperature one can read off the 
intermonopole potential from eq. (|52|) . The potential is the Debye screened 3D Coulomb 
potential. At nonzero temperature eq. (B3) gives it as the 2D Coulomb potential at short 



8 Our calculation is only valid in the range of temperatures below the Higgs expectation value h. At 
T > h the theory becomes essentially nonabelian, since the Higgs particle as well as W bosons are present 
in the "vacuum" in large numbers, and our methods become inapplicable. In fact in this range there may 
be another phase transition corresponding to the "melting" of the Higgs expectation value. Our analysis 
in this and previous sections has nothing to say about this regime. 

9 Strictly speaking V creates a monopole - that is the monopole with half the magnetic charge than 
the 'tHooft-Polyakov solution. However the interaction between these objects is qualitatively the same. 
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distances, again screened at distances x 2 > At high temperature then we see from eq. 
(|54j) that the interaction between the monopoles is linear with the "string tension" equal 
to M s . Thus the behaviour of the monopoles is in many senses "dual" to that of electric 
charges^- 

Finally we wish to comment on the behaviour of the spatial Wilson loop in this 
model. One expects the spatial Wilson loop to have an area law. This area law has a 



simple interpretation in terms of the thermal ensemble of magnetic vortices |T7||. The 
vortex picture is generally true and as discussed in |J, holds even in theories where the 
Wilson loop at zero temperature has a perimeter behaviour. We have nothing to add to 
the discussion of the vortex ensemble picture, since our present approach is not formulated 
directly in these terms. It is curious to note however, that the spatial string tension can not 
be calculated at all in the dimensionally reduced theory of the field x-, eq. (|27|) , although it 
can be easily shown that it does not vanish. Consider the calculation of the Abelian Wilson 
loop W(C). As discussed in ||, it is the generator of the magnetic Z 2 transformation inside 
the contour C. The path integral representation of this expectation value is 

(d^x ~ U) ( d »X ~ jp) + + cos 2 X \ (57) 



< W (C) >= jd X ex V | - J d 3 x£^ 



with 

j^x) = tt5^5(x ), xeS (58) 
= 0, x <£ S (59) 

and S is the area bounded by C lying in the xq = plane. This expression clearly can not 
be dimensionally reduced, since the external current has components at all Matzubara 
frequencies. If we nevertheless naively disregard the higher Matzubara modes of the x 
field, we find|^| 

W{C) = e~£ s (60) 



10 Of course one always has to keep in mind that physically monopoles and charges are very different 
objects in this models: the charges are particles, while the monopoles are instantons. 

11 The contribution of the zero Matzubara frequency is trivial, since the path integral in this sector is 
dominated by the configuration 4>{x) — 7r, x £ S which has vanishing action. 
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Since the cutoff on the reduced theory, a is of order 1/T, parametrically this result is 
precisely what we expect a oc g 2 T. However to determine the exact numerical value we 
would have to look into details of the nonzero Matsubara mode configurations. 

6 Conclusions 

In this paper we have studied the deconfining phase transition in the 2+1 dimensional 
Georgi Glashow model. This model to this day remains the only analytically solvable 
confining theory, and it is therefore interesting to see how the deconfining transition occurs 
in this simple yet nontrivial setting. 

We found that the phase transition is indeed associated with the restoration of the 
magnetic Z 2 symmetry with the universality class of 2D Ising model. We have shown 
that the transition is driven by the appearance of the charged plasma rather than by the 
binding of the monopoles. The presence of the monopoles is significant, since it pushes the 
value of the critical temperature up by a factor of two (in the infinitely heavy Higgs limit) 
relative to the plasma transition in the noncompact electrodynamics. We have derived the 
effective action in the high temperature phase without recourse to the derivative expansion 
and have calculated the correlation function of the vortex operators. These correlators 
decay exponentially above the phase transition, which is consistent with the restoration of 
the magnetic symmetry. This also shows that the monopoles are bound linearly at high 
temperatures due the charged plasma effects, and these effects are much more important 
than the monopole self interaction, which by itself leads only to a logarithmic intermonopole 
potential. 

We believe that these qualitative features are universal and should generalize not 
only to the strongly interacting non Abelian theories in 2+1 dimensions, but also in large 
measure to 3+1 dimensional QCD along the lines discussed in ||. Unfortunately we are not 
yet able to generalize the present quantitative analysis to these more interesting theories. 
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